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Galaxy cluster outskirts: a universal entropy profile for relaxed 
clusters? 
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ABSTRACT 

We fit a functional form for a universal ICM entropy profile to the scaled entropy profiles of 
a catalogue of X-ray galaxy cluster outskirts results, which are all relaxed cool core clusters 
at redshift below 0.25. We also investigate the functional form suggested by Lapi et al. and 
Cavaliere et al. for the behaviour of the entropy profile in the outskirts and find it to fit the 
data well outside 0.3r 2 oo- We highlight the discrepancy in the entropy profile behaviour in the 
outskirts between observations and the numerical simulations of Burns et al., and show that 
the entropy profile flattening due to gas clumping calculated by Nagai & Lau is insufficient 
to match observations, suggesting that gas clumping alone cannot be responsible for all of the 
entropy profile flattening in the cluster outskirts. The entropy profiles found with Suzaku are 
found to be consistent with ROSAT, XMM-Newton and Planck results. 
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1 INTRODUCTION 

A breakthrough in the study of the low surface brightness X-ray 
emission from the intracluster medium (ICM) in the outskirts of 
galaxy clusters has been made possible by the low and stable parti- 
cle background of Suzaku, and a statistical sample of clusters stud- 
ied to has started to emerge (see table Q}. All of these are 
relaxed cool core clusters with redshifts of less than 0.25. Since 
these clusters have been chosen because they are very X-ray bright, 
they do not provide a representative sample of the actual cluster 
population, but nonetheless general trends have started to emerge. 
Understanding the thermodynamic properties of clusters in the out- 
skirts is important for constraining the physical processes occur- 
ring in the ICM. We also need to understand the radius at which the 
hydrostatic equilibrium approximation breaks down, and the con- 
tributions of non-therma l pressure suppo rt which are expected to 
increase in the outskirts dLau et al]|2009l) in order to calculate ac- 
curate masses out to r2oo and beyond. 

One of the interesting findings from studies of the outskirts of 
galaxy clusters is that the scaled ICM entropy (S = kT/n^/ 3 ) 
profile has a similar shape for each cluster, flattening off above 
~ 0.5r2oo away from a powerlaw increase. The entropy profile is 
of fundamental importance because it both determines the structure 
of the ICM and provides a record of the t hermodynamic history of 
the ICM. As described in lVoit et"al]d2005l) . the entropy is more rep- 
resentative of the thermodynamic history of the ICM than temper- 
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1 7"200 is the radius within which the mean mass density is 200 times the 
critical density for a flat universe. 



ature because when heated gas expands in a gravitational potential 
its thermal energy can be converted into gravitational potential en- 
ergy. Introducing heat will always cause the entropy to rise while 
losing heat through radiative cooling will always cause the entropy 
to fall, whereas it is possible for the luminosity-weighted tempera- 
ture to rise only slightly on the input of a large amount of ener gy, 
making it a less sensitive probe of ICM history dVoit et ail2002l) . 

Here we compile a sample based on the papers summarised in 
table [T] and fit an analytic function to the scaled entropy profiles. 
In addition to the clusters studied with Suzaku, we also add Abell 
1835 and Abell 2204 which have been studied with Chandra out to 
rioo and 0.7r2oo respectively. We also add the Virgo cluster which 
has been studied with XMM-Newton out to r^oo along one nar- 
row strip. We do not i nclude the clusters studied with Suzaku in the 
outskirts described in Aka matsu & Kawaharal d201ll) as these ob- 
servations are aligned specifically with radio relics, and so are not 
representative of the ICM of the whole cluster. Compiling a sam- 
ple like this reduces the effect of the small azimuthal coverage of 
some of the observations, allowing a more realistic average to be 
obtained which reduces the bias of looking along particular direc- 
tions in a cluster w hich may not be repr esentative of the cluster as a 
whole. For instance lEckert et alj d2012h has used ROSAT PSPC ob- 
servations and found that clusters in the outskirts have significant 
(~ 70 percent) azimuthal scatter in their surface brightness when 
the clusters are divided into 12 sectors. Theoretically galaxy clus- 
ters are expected to be increasingly anisotropic in their outskirts 
with accretion o ccurring preferent ially along large scale structure 
(LSS) filaments dBurns et alfcOld) . 

As mentioned in Walk er et alj d2012h b, the emerging picture 
for these low redshift, relaxed cool core clusters is that of an en- 
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Table 1. Sample of galaxy cluster outskirts observations used. 



Cluster z Reference Plot symbol 



which can be written as: 

SyS(0.3r 200 ) = A c (r/R) Bc e Cc(1 - (r/H)) (1) 

where R is the virial radius, taken to be half the turnaround ra- 
dius for the collapsing overdensity (which has detached from the 
Hubble flow) from which the cluster forms, which in a ACDM 
universe is R — rioo- This functional form has a m aximum at 
r/R = B g /Cc, and a gradient a t R which is Be — Cr-. lLapi et al.l 
(2010) and lCavaliere et aljfeoilh propose that as clusters get older 
the flattening and downturn of the entropy profile at large radius 
becomes more severe due to the reduction in the strength of the 
accretion shock. However the small redshift range covered by our 
sample (all have redshifts less than 0.25) means that no age evolu- 
tion is discernible within our sample. The entropy profile flattening 
is not expected to occur in non-cool-core clusters, since for these 
younger clusters the high accretion rates and strong shocks should 
allow the entropy profile to follow the r 1 ' 1 relation out to the virial 
radius, however their often disturbed morphologies complicates ob- 
servations. 



2 ANALYSIS 

We use the entropy profiles for the clusters listed in table [TJ and 
scale them by their entropy at 0.3r2oo — ?"2500, a region where the 
clusters all demonstrate a powerlaw increase in entropy and are in 
hydrostatic equilibrium since this is both far enough away from the 
core to avoid the affects of cooling and feedback processes, and 
far enough from the outskirts to be unaffected by accretion pro- 
cesses . This is shown in Fig.[TJleft. For the data from lBautz & et all 
(2009), which presents analytic forms for the entropy profile in- 
stead of points, we show points centred on the annuli from which 
the spectra were extracted, and include only the data in the range 
in which the analytic forms for the temperature and density (and 
thus entrop y) were fit (between 5 '-ll.5'). We do the same for the 
results of iBonamente et al.l ( 120 12I) . in which a best fit analytic en- 
tropy profile line is presented rather than entropy data points for 
each annulus. 

The results for the Virgo cluster from lUrban et all ( l201ll) . 
shown by the green crosses in Figs.Q]and[3] are systematically be- 
low the general trend of the other clusters, indicating more severe 
entropy profile flattening, with significant oscillatory behaviour in 
the entropy profile with radius. This may be the result of a combina- 
tion of the following factors. Virgo is the poorest and least massive 
cluster in our sample, which may indicate that entropy profile flat- 
tening is a function of cluster mass. The small azimuthal coverage 
of the Virgo observations (as shown in Fig. [2] at 7-200 the azimuthal 
coverage is ~ 2 percent) may mean that the results are not repre- 
sentative of the cluster as a whole, given the observed azimuthal 
variatio ns in the surface br ightness of the ICM in the outskirts of 
clusters teckert et alj2012t) . For instance, if the observed strip hap- 
pened to coincide with a region of higher than average gas clump- 
ing then this would anomalously give a larger inferred density, and 
thus a lower inferred entropy than the actual azimuthal average in 
the outskirts. The fact that there is a cold front at ~ 0.2r2oo in the 
direction observed may also make the results unrepresentative of 
the whole cluster. 

The outermost entropy measurement from Abell 1413 (blue 
triangles) also lies above the average relation. However as shown 
in Fig.|2]the fraction of the ICM in the annulus which is studied is 
quite low (~ 15 percent), and so may not be representative of the 
whole cluster at that radius. 
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tropy profile which obeys the r 1 ' 1 powerlaw (from lvbit et alj2005T) 
from around 0.2r2oo out to around 0.5r2oo, beyond which it flat- 
tens. The radius at which this flattening starts, and the extent of this 
flattening, provides an important test of numerical simulations of 
galaxy clusters, requiring an accurate understanding and modelling 
of the gravitational and non-gravitational heating mechanisms in 
the ICM. 

A number of physical processes are expected to affect the en- 
tropy profile in the outskirts. Gas clumping causes the gas density 
to be overestimated, thus causing the entropy to be underestimated. 
Gas c lumping however is expected to be most significant outside 
r 20 o, JNagai &Laull201ll) . 

A difference between the electron and ion temperatures in 
the ICM inside the accretion shock i n the outskirts may als o be 
important, and has bee n proposed in Hos hino & et al.l d2010h and 
I Akamatsu et al.l d20Hh as a cause for entropy profile flattening. 
This would arise because the ions, which carry most of the kinetic 
energy of the accreting gas (due to their greater mass) are ther- 
malised first after the accretion shock, and this energy is then trans- 
ferred inefficiently and slowly to the electrons through electron-ion 
collisions. However the absence of a pressure drop in the outskirts 
of galaxy clusters in the Sunyaev-Zeldovich (SZ ) effect observa- 
tions with Planck ( Planck Collab oration et al .120121) , and the agree- 
ment between the Planck observations of the gas pressure with 
simulations suggests that the electrons and ions in the outskirts 
are in equilibrium. In a ddition, the hydrodynamical simulations of 
Wong & Sarazi r\ d2009l) indicate that within rioo the electron and 
ion t emperature s are e xpect ed to differ by less tha n a percent. 

lLapi et all |2010) and ICavaliere et al] d201 ll) have proposed 
that entropy profile flattening in the outskirts of relaxed galaxy clus- 
ters is the natural result of the weakening of the accretion shock as 
it expands. This is a result of the accreting material originating from 
the tapering wings of the dark matter perturbation over a decreasing 
background density level due to the accelerated cosmic expansion. 
Therefore the accreting gas falls through a progressively smaller 
potential drop as the accretion shock expands outwards. The weak- 
ening of the accretion shock as it expands reduces the entropy gain 
at the shock and also increases the amount of bulk energy passing 
across the shock by the infalling gas. This increases the level of tur- 
bulence and non-thermal pressure support in the outskirts, causing 
hydrostatic equilibrium mass estimates which use only the ther- 
mal pressure to underestimate the true mass (as has been found in 
lLau et alj|2009h . 

Lapi et al.l (2010) and ICavaliere et alj d201 ll) propose a func- 
tional form for the entropy profile in the outskirts (outside ~ R/3) 
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Figure 1. Le/hEntropy profiles for the clusters sh own in tab le Q] scaled by S'(0.3r20o) ■ Individual clusters are colour coded as shown in table[T] The solid 
black line shows the r 11 powerlaw relation from lVoit et alj 120051) . Righf.We plot S(r)/r (scaled to 0.3r2on) to show the deviation from a powerlaw more 
clearly. The black line is the best fit lin e to the data outside 0.2r200 using a form S'/S(0.3r20o) = A(r/r20o) 1 ' 1 e~ ^ '/ Br wo) The best fit using the 
functional form of Cavaliere e t al.l 1201 ll) (equation[TJ is shown by the blue line. For each model the 2 a variations calculated using Monte Carlo methods are 
shown by the dash ed lines. The solid red lines show the range produced by density variations of 30 percent, which is the observed azimuthal density variation 
found near T200 m Eckert et al, 12 0121) . 
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Figure 2. Percentage azimuthal coverage as a function of radius for the 
observations used. 



In Fig. [TJ (right) we plot S/r against r (scaling the pro- 
files by 5'(0.3r2oo)/0.3r2oo), which more clearly shows the de- 
viation from a simple powerlaw above 0.5r2oo. We find that the 
profile is fitted well by the functional form S'/S(0.3r2oo) = 
^(r/r2oo) 11 e~ (r/Sr20o)2 for r ^ 0.2r 2oo with best fitting pa- 
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We also find the best fit to the scaled entropy profiles in the 
range r ^ 0.3r2 oo us ing the function a l form of equation [TJ from 
lLapi et al.l j201Ch and ICavaliere et all ( 1201 ll) . which is found to 



model the entropy profiles well with best fit parameters Ac 



1.02Z 



B c = 1.8t u LC c = 3.3 



+ Q 2, so the best fit relation is 



S/S(0.3r 20 o) = 1.02(r/i?) 1 ' 8 e 3 ' 3(1 - (r/fl)) 



(3) 



Since the errors on each parameter are correlated, the errors 
on the best fits were obtained by using a Monte Carlo method with 
10000 trials, and the 2 a variations of the best fit models are shown 
by the dashed lines in Fig. [T] right. Black lines show equation [2] 
while the blue lines show equation [3] When performing the fitting 
the entropy profiles from each cluster were also weighted by the 
azimuthal coverage of the observations of each cluster (shown in 
Fig. [2}, so that more weight was given to observations with larger, 
more representative azimuthal coverage. This reduces the possible 
bias of observations which were taken along narrow strips which 
may not be representative of the cluster as a whole. 

The solid red lines in Fig. [T] (right) show the effect of 30 per- 
cent density variations on the best fit entropy profile. This is the 
level of azimuthal scatter in the gas density inferred from the az- 
imuthal scat t er in t he surface brightness of the clusters studied in 
lEckert et alj 1201 2h (where the observed surface brightness scatter 
was ~ 70 percent around r2oo). We find that the majority of the 
data lie within this range around the best fit profile, suggesting that 
most of the scatter around the best fit profile can be explained by 
the ~3 percent azimuthal density variations found in Eck ert et"al] 
( 120121) . The Virgo results are however inconsistent with the trend of 
the othe r clusters. Th i s may be because the azimuthal scatter mea- 
sured in lEckert et alj d2012l) was found by dividing the clusters in 
their ROSAT sample into 12 sectors of opening angle 30 degrees, 
whereas the Virgo strip is much narrower than this (its opening an- 
gl e is ~ 8 degrees). It is therefore possible that the scatter measured 
in lEckert et alj J2012h underestimates the level of scatter at scales 
smaller than the sector size they used. 

In Fig. [3] (black lines) we compare the scaled entropy profiles 
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Figure 3. Same as F ig. [T] left but including the results of the simulations 
of Burns et al. ( 2010) as the black solid line. The dashed black lines show 
th e standard deviat ion from the sample of numerical simulations studied 
in Burns et alj <2010l) (which were normalised at 0.5r2oo for the sample, 
hence the zero scatter around t his point). The blue solid line shows the re- 
sults of the CSF simulations of Nagai & Lau 1 201 1) showing the flattening 
of the entropy profile as a result of gas clumping. In both cases the level of 
entropy flattening and turnover, and the radius at which the entropy starts to 
turnover, disagree with the observations. 



with the similarly sca led predicted ent ropy profiles from the nu- 
merical simulations of lBurns et al.l bold) obtained using the Eule- 
rian adaptive mesh refinement (AMR) cosmology code Enzo using 
only adiabatic gas physics (thus neglecting non-gravitational heat- 
ing, cooling, magnetic fields, and cosmic rays). We find that the ra- 
di us at which entropy flattening begins is higher in the simulations 
of lBurns et all l feOlOl) than in the observations, suggesting that more 
physical processes need to be included in the numerical simulations 
to accurate ly describe the observed pro perties of clusters. The sim- 
ulations of iKravtsov & Borganil d2012l) also show the entropy pro- 
file flattening starting outside rino, at a higher radius than is seen in 
the observations. A possibility is that the simulations may overesti- 
mate the number of major mergers at late cosmic times, causing the 
strength of the accretion shocks to be overestimated and increasing 
the radius at which the entropy profile flattening begins. 

In Fig. [3] we also compare the observed entropy profiles with 
the profile shape o btained in the gas clumping simulations of 
iNagai & Laul feoill) (solid blue line), which included the effects 
of cooling and star formation (CSF). The clumping effect only be- 
comes significant outside r2oo, and is not large enough to account 
for the amount of entropy flattening and turnover. 



3 COMPARING TO ROSAT, XMM-NEWTON AND 
PLANCK RESULTS 

IPlanck Collaboration et al. (2012) obtained a pressure profile for 
a stacked sample of 62 clusters using the S Z-effect reaching out 
to 3r5oo — 2r2oo, while lEckert et al.l ( 12012) obtained an average 
density profile from a sample of 31 clusters obtained with ROSAT 
PSPC observations out to r2oo- Using that S oc P/ri^J 3 we can 
combine these to obtain a scaled entropy profile, and this is plotted 
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Figure 4. 7bp:The thick solid black line shows the entropy prof i le ob tained 
by combining the pressure profile of Planck Collaboration et al. (20l2) (ob- 
tained from a sa mple of 62 clusters) and the density profile obtained with 
ROSAT PSPC in lEckertetaTI <2012h (obtained from a sample of 31 clus- 
ters). The dashed thick black lines show the scatter. Bottom: The solid blue 
line shows t he entropy p r ofile o btained by combining the ROSAT density 
profile from Eckert et al. (2012) with the extrapolated XMM-Newton tem- 
perature profile shown in Fig. [5] and the \-a error is shown by the thick 
blue dashed lines. In both plots the thin black line shows the r powerlaw 
relation. 



in Fig. [4] (top) as the thick black line. The resulting entropy pro- 
file agrees well with the Suzaku observations, flattening off above 
0.5r 2 oo- 

As another test we use t he average temperature profile found 
in lLeccardi & Molendil j2008l) out to ~0.7r2oo using XMM-Newton 
observations of a sample of ~ 50 clusters, and extrapolate this 
out linearly to noo, as sho wn in Fig. HI When combined with the 
ROSAT density profile from Eck ert et alj J2012h the resulting scaled 
entropy profile is shown by the solid blue line in Fig. [4] (bottom), 
which again agrees well with the Suzaku results. 

In order for the entropy profile to follow a r relation out to 
r2oo and beyond, the temperature profile would need t o be the blue 
line shown in Fig. \5\ when the density profile from ( Eck ert et alj 
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I2012I) is used, which shows a much shallower temperature decline 
in the outskirts than is observed. 



4 SUMMARY 

We have compiled a sample of relaxed cool core galaxy clusters 
(with z ^ 0.25) whose ICM has been studied to ~ r2oo- We 
find that the scaled entropy profiles are best fit using a function of 
the form S/S{0.3r 20 o) = 4.4(r/r 2 oo) 11 e ~ (r/r20o)2 . We also iri - 
vestigate the analytic function proposed bv lCavaliere et al.l ( l201lh . 
and have found it to fit the data well out side 0.3r2o o (eq uation 
3). This adds support to the suggestion of lLapi et alj d2010h and 
Cavali ere et al. \ 1) that the flattening and downturn of the ICM 
entropy in the outskirts is the result of the decreasing strength of the 
accretion shock as it expands, due to the accreting matter coming 
from the tapering wings of the dark matter perturbation. The in- 
crease in bulk energy seeping through the shock is expected to lead 
to turbulence and non-thermal pressure support, which are predic- 
tions that could be tested by next generation instruments studying 
cluster outskirts. 

The observed scaled entropy profiles turnover at a radius lower 
than t hat predicted by the numerical simulations of Bur ns et al.l 
(2010) which used only adiabatic gas physics. The entropy profile 
flattening d ue to gas clumping calculated in the CSF numerical sim- 
ulations of iNagai & Laul d201ll) is insufficient to match the obser- 
vations (Fig. [3}, suggesting that clumping alone cannot be respon- 
sible for the flattening and turnover of the entropy profile, indicat- 
ing that some other process must contribute. Clumping is required 
to explain the observed excess of the gas mass fraction over the 
mean cosmic baryon fraction around r?oo {simionescu et al J20T1I ; 
IWalker et al.ll2012b : IFlanck Collaboration et alj2012h . though part 
of this excess may be the result of increasing non-thermal pressure 
support in the outskirts causing the total mass estimates using only 
thermal pressure to be underestimat ed by around 20 percent at r2oo 
dLau et all2009t IWalkeret alj2012b ). 

The Suzaku entropy profiles are in strong agreement with the 
entropy profiles obtai ned by combining the mean press ure profile 
obtained with Planck (Planck Collaboration et al. 2012), the mean 
density profile obtained with ROSAT PSPC feckert et al.l2012h and 
the extrapolati on of the mean temperatu re profile obtained with 
XMM -Newton fceccardi & Molendill2008l) . 
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